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We use the numerical optimization techniques of Uskov et al. [l[ to investigate the behavior of 
the success rates for KLM style [2|] two- and three-qubit entangling gates. The methods are first 
demonstrated at perfect fidelity, and then extended to imperfect gates. We find that as the perfect 
fidelity condition is relaxed, the maximum attainable success rates increase in a predictable fashion 
depending on the size of the system, and we compare that rate of increase for several gates. 
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PACS numbers: 42. 50. Ex, 03.67.Lx, 42.50.Dv 



I. INTRODUCTION 

Optics provides one of the most promising implemen- 
tations for quantum information processing, due to long 
photon decoherence times, the ease of photon manipu- 
lation, and the ability to transmit quantum state in- 
formation over very large distances. Optics have also 
been suggested as buses in hybrid matter-optical quan- 
tum computers. It is therefore desirable to be able to effi- 
ciently create and manipulate states at the single photon 
level. In this work, we will focus on manipulating states 
through optical gate implementation, rather than on sin- 
gle photon or quantum state generation |3[. 

Any number of implementation schemes have been pro- 
posed to preform non-trivial operations at the single pho- 
ton level. These include non-linear materials, coupling 
through atomic interactions, Zeno-type non-unitary in- 
teractions |4[ , and silica-on-silicon waveguides [5[ . 

Knill, Laflamme, and Milburn significantly advanced 
the prospect of single photon quantum computing in their 
seminal paper [2], in which they overcame the need for 
nonlinear interactions by using the inherent nonlinearity 
of photon measurements. In this scheme, the compu- 
tational system is combined with ancillary modes, and 
the gate operation is performed on the enlarged state 
space, red states of the computational system when 
the measurement is successful). The ancilla modes are 
measured with photon- number-resolving detectors, leav- 
ing the computational modes undisturbed and in the de- 
sired output state provided the measurement is success- 
ful. The probabilistic nature of quantum measurement 
implies a trade-off between the success rate of the oper- 
ation (the probability of obtaining the desired measure- 
ment outcome for the ancillary modes) and the fidelity 
(the overlap between the actual and desi 

In previous work [1], the authors have shown that 
a combination of analytical and numerical techniques 
may be used to design optimal linear optics transfor- 
mations implementing two- and three-qubit entangling 
gates. Specifically, solutions were obtained that have 
the maximum possible measurement success probability 



given the constraint of perfect fidelity for a given desired 
gate. In practical implementations, however, the goal 
of perfect fidelity may not always be desirable or even 
attainable (for given ancilla resources). In the present 
work, we therefore generalize our previous techniques to 
the case of imperfect fidelity, and investigate the above- 
mentioned trade-off between the fidelity and success of 
the linear optical transformation. We will show that for 
sufficiently small deviation from perfect fidelity, a single 
parameter determines the relationship between fidelity 
and optimal success rate. We will also observe that once 
the perfect fidelity condition is relaxed, the resulting rel- 
ative gain in success rate appears to grow with the size 
of the linear optical system under consideration. 



II. LOQC GATE OPTIMIZATION 

A Linear Optical Quantum Computing (LOQC) 
measurement-assisted transformation, as suggested by 
KLM, is schematically illustrated in Fig. [TJ Here the in- 
put state |^ mp ) x |^ ancilla ) is a product of the computa- 
tional state and an ancilla state, with |\^° mp ) containing 
M c photons in N c modes, and |\[/ ancilla ^ containing M a 
photons in N a modes. 

The N c computational modes are those on which the 
actual gate is intended to act. Assuming the Fock basis 
and dual-rail encoding, each qubit is represented by a 
single photon in two computational modes: specifically 
the logical states 1 1) and \\) are represented by | li, O2) 
and |0i, I2), respectively. The number of computational 
qubits is then equal to the number of photons M c in the 
input state, which is half the number of modes: M c = 
N c /2. This also implies that any output state with more 
or fewer than one photon in any logical qubit is not a 
computationally valid output. 

In Sec. lIII[ we first discuss two-qubit gates, so the com- 
putational state will consist of M c = 2 photons in N c = 4 
modes. However, our numerical optimization approach 
applies to any gate operation [1], e.g., the three-qubit 
Toffoli gate investigated in Sec. IIV( where the computa- 
tional state consists of M c = 3 photons in N c = 6 modes, 
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and others. 

The ancilla state may be chosen from a wide array 
of possible states, and this state may be separable, en- 
tangled, or an ebit state carrying spatially distributed 
entanglement [gj]. We note that there is a lower limit 
to the amount of ancilla resources needed to implement 
any given gate. Additionally, we allow for the possibil- 
ity of an arbitrary number N v of vacuum modes that are 
unoccupied in the input state. 
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FIG. 1: (Color Online) A scheme for an LOQC transformation 
or a quantum state generator pj. The computational input 
state is a separable state of two or more dual-rail encoded 
qubits. The ancilla state may be separable or entangled. 

The linear optical device illustrated in Fig. [I] trans- 
forms a set of input creation operators a\ m ^ to a weighted 

sum of output creation operators J2j Uija^ ut ^ . Here U, 
which contains all physical properties of the device, is 
an N x N unitary matrix, where N = N c + N a + N v is 
the total number of modes (computational, ancilla, and 
vacuum). The total input state may be written in the 
Fock representation as |^i n ) = ^2, . . . , wjv), where rti 
is the occupation number of the i-th input mode, and 
Tii = M c + M a = M is the total number of photons. 
The transformation Cl then takes the form 



|*out> = n|*i» 
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We note that Cl is a multivariate polynomial of degree M 
in the elements Uij. 

Once the transformation is complete, a measure- 
ment is applied to the N a + N v ancilla and vacuum 
modes. This measurement can be formally described 
by a Kraus POVM Q operator acting on these modes 

only: P = \0n c +i, 0iV c +2, On) (^measured |. In the 
case of a number-resolving photon-counting measure- 
ment, (^ m easured| = (&jv c +i , k Nc + 2 , • • • , k N | , where ki 
is the number of photons measured in the i-th mode 
of the ancilla. Finally, the resulting transformation of 
the computational state is a contraction quantum map 

i*~r> = A\^r p )/\\Amr p )\\ m, a = aw is 

defined by 



The linear operator A contains all the information of rel- 
evance to the gate or state transformation. Again the 
entries of the A matrix, Aij = Aij(U) are multivariate 
polynomials of degree M in the entries of U . 

We note that the dual-rail computational basis is a sub- 
set of all possible states of M c photons in the 2M C com- 
putational modes. Since the input state of the compu- 
tational modes necessarily belongs to the computational 
space, while the output state may in general be any state 
of M c photons in 2M C modes, the transformation matrix 
A is a rectangular matrix, mapping the computational 
Hilbert space, of dimension 2 Mc , to a larger Hilbert space, 
of dimension (3M C - 1)!/(2M C - 1)!(M C )!. For example, 
A is a 10 x 4 matrix for a two-qubit gate. Similarly A 
is a 56 x 8 matrix for a three-qubit gate and 35 x 9 for 
two-mode biphotonic qutrit gates [8]. 

We now define precisely the operational fidelity of a 
transformation, which in general differs from the com- 
mon measure of fidelity for a state transformation [gj]. 
Physically, the transformation A has 100% fidelity if it 
is proportional to the target transformation A Tar , i.e., 
A = gA Tajr , where g is an arbitrary complex number 
(in which case S = \g\ 2 is the success probability of the 
transformation [lQ[). In general, we may consider com- 
plex rays f3\A and /32A Tar , /?i,/?2 C C, as elements of a 
complex projective space, and define the fidelity as 



F(A) 



(i|i Tar )(i Tar |i) 
<A|A)(i Tar |i Tttr ) ' 



(3) 



where the Hermitian inner product is (A\B) = 
Tr(^ B)/D c , and D c = 2 Mc is the dimension of the com- 
putational space. F is closely related to the Fubini-Study 
distance 7 = cos~ 1 (^/F) [ll[, but for numerical compu- 
tations F has the advantage of being non-singular near 
F = 1. 

In general, the success probability S depends on 
the initial state |\I/^ mp ). S is bounded above by the 
square of the operator norm, ||A|| 2 = (||A|| Max ) 2 = 
Max((^^ n omp |iti|^ n om P)) 5 an d below by (||i|| min ) 2 = 
Min((^° mp |iti|^° m P)) j where the maximum and min- 
imum are taken over the set of properly normalized in- 
put states. As a more convenient measure, we use the 
Hilbert- Schmidt norm It is easy to verify that 

(||i|| Min ) 2 < (A\A) < (||i|| Max ) 2 . As fidelity F -+ 1, 
||i|| Min /||i|| Max -> 1 and S becomes a well-defined quan- 
tity equal to the Hilbert-Schmidt norm. In the following, 
we refer to 



S(A) = (A\A) 



(4) 
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as the success probability, keeping in mind that it may 
not correspond to the success probability for every initial 
state, except in the case of perfect fidelity. Since the 
state transformation A is a function of the linear mode 
transformation U, the fidelity F and success S are also 
functions of U. 
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In the above discussion, we have assumed unitarity 
of the underlying linear mode transformation U. How- 
ever, an arbitrary (N c -\-N a ) x (N c -\-N a ) complex matrix 
U of unit spectral norm (||?7|| = 1, where \\U\\ is the 
largest singular value of U) may be shown via the uni- 
tary dilation technique [HI, G3 to be equivalent to an 
(N c + N a + N v ) x (N c + N a + N v ) unitary matrix having 
the same success and fidelity, where N v is the number 
of vacuum modes in the input and output states. Thus, 
it is very convenient in practice to relax the unitarity 
condition and consider general matrices U of norm 1, 
with the understanding that the number of singular val- 
ues in U different from unity corresponds to the number 
of vacuum modes that are required to implement that 
solution [l3|| . Furthermore, we may consider arbitrary 
complex matrices U with the rescaling U —> U/\\U\\ to 
ensure unit norm. The fidelity function (j3j) is unaffected 
by the rescaling, while the success function (jl]) must be 
generalized as 

S(A) = (i|i)/||C/|| 2M (5) 

for general complex U. The success S(A) exhibits cusp- 
like singular behavior in this extended search space when- 
ever the largest singular value of U goes through a double 
or higher-order degeneracy. This occurs in particular in 
the neighborhood of the manifold WU = I. 

The optimization problem we address is to maximize 
the success probability S for a given target transforma- 
tion A Tar , for given ancilla resources, and for a given 
fidelity level F < 1. This is mathematically equivalent 
to unconstrained maximization of the function S + F/e 
in the space of all N x TV matrices U, where 1/e is a La- 
grange multiplier. Here e — >• + corresponds to maximiz- 
ing the success probability while requiring perfect fidelity 
(F = 1). As e is increased, the maximum of S+F/e yields 
linear optics transformations that maximize the success 
S as a function of the fidelity F. 

Given one optimal transformation U that (locally or 
globally) maximizes success S for a given fidelity F, e 
may in general be continuously varied to obtain a one- 
parameter family of optimal transformations, tracing out 
a curve in success-fidelity space. Multiple such families 
of solutions may coexist, corresponding to different local 
maxima of the success rate at a single given value of the 
fidelity. 

Numerically, as we attempt to follow a family of op- 
timal solutions, we find that it is necessary to perturb 
the U matrix slightly at each step before incrementing e 
and repeating the optimization process. This is due to 
the fact that such a family is typically associated with a 
degeneracy of several, or in some cases all, singular val- 
ues of U. The family of optimal solutions then exists 
along a multi-dimensional ridge of the success function 
S (Eq. (j5j)), with cusplike behavior of S and thus of the 
optimization function S + F/e along directions in U space 
orthogonal to the ridge orientation. A slight perturba- 
tion of the U matrix at each successive value of e moves 



the solution away from the cusp and allows the ridge to 
be followed in a continuous manner. 



III. TWO-QUBIT GATES 

In this section we examine two controlled-phase gates, 
CNOT (0 = 180) and cs(90°) (0 = 90), which are 
of course equivalent to more general controlled-unitary 
gates by local qubit transformations. We also consider 
the more general B gate, which is notable because an ar- 
bitrary two-qubit SU (4) gate may be constructed out of 
two B gates and local rotations, as shown by Zhang et 
al. [14[. The maximum success rates for the two phase 
gates was first determined numerically by E. Knill [12|. 
Subsequently the maximum success rates of all of the 
controlled-phase gates and most general gates were found 
numerically by the authors These solutions were 

found at perfect fidelity. A question that remains to be 
answered is: what are the maximum success rates for 
imperfect gates with fidelity (slightly) less than 1? 

We use the Cartan decomposition KAK to define each 
of the target transformations (cs(90°), CNOT, and the B 
gate) using only three real parameters (c\. C2, C3). For 
details on the Cartan decomposition see [20 and 
In this decomposition, CNOT (also equivalent to CZ via 
local rotations) is represented by the point (7r/2,0,0), 
cs(90°) is represented by (7r/4,0,0), and the B gate is 
represented by(7r/2, 7r/4, 0). 



Success Rate of the CNOT gate vs. 8 
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FIG. 2: (Color Online) Numerical results for the success rate 
S of the CNOT gate as a function of S = 1 — F for our per- 
turbative relaxation of the F — 1 constraint. The continuous 
curve is a fit given by Eq. ©. 

We begin with the CNOT gate as shown in Fig. O We 
note that M c = 2 computational photons in N c = 4 com- 
putational modes are required for any two-qubit gate in 
the dual rail representation. Additionally, the minimum 
ancillary resource required to implement CNOT with per- 
fect fidelity consists of two single-photon ancillas, i.e. the 
CNOT gate requires M a = 2 ancilla photons in N a = 2 
ancilla modes [13| . Therefore we have a total of four pho- 
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tons in six modes in the initial state, and the linear opti- 
cal transformation U is a 6 x 6 matrix. Since the fidelity 
F is a smooth function of U, the fidelity behaves quadrat- 
ically, F = 1 — 0(r] 2 ), around the maximum F = 1, while 
the success S is in general not at a local maximum at 
the F = 1 point and therefore behaves generically as 
5 = So + ar] + fr?? 2 + • • • . Here 77 is a local coordinate 
in U space along the direction of the gradient of S (with 
r] = corresponding to a solution with perfect fidelity), 
So is the success rate for perfect fidelity, and a — |VS|. 
We then immediately see that the curve of success vs. 
fidelity should be of the form, 



S(F) = So + Si(l - F) 1/2 + S 2 (l - F) 



(6) 



in the area around F = 1, i.e., a power series expansion 
in (1 — F) 1 / 2 . Defining 5 = 1 — F, we have 



S(5) = S + S 1 S 1 / 2 + S 2 S- 



(7) 



Looking at Fig. [2j we see that the data from our numer- 
ical optimization does follow this form very well, with a 
best fit 



S(S) = 0.0740 + 0.0765 5 1/2 + 0.0199 S . 



(8) 



the Hilbert-Schmidt norm, i.e., these are the success rates 
averaged over all possible computational input states. 
We find that while the spread between the minimum and 
maximum success rates is < 10 -4 for the 3rd family of 
solutions, the corresponding spread for the Knill family 
and 2nd family of solutions is of the same order as the 
average success rate. This suggests that in applications 
where single quantum state generation or transformation 
of a single input state is desired, the method described 
here may be used to design a linear optical device that 
maximizes the success rate for that specific application. 

Several Families of Transformations for the CNOT gate vs. 6 
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When F = 1, S = and the function S(S = 0) reduces 
to the S = 2/27 « 0.0740 result found by Knill [H, 
which was confirmed previously by the authors [1]. The 
ratio Si/ So = 1.03 contains the most interesting informa- 
tion about the system, as it is a measure of the relative 
increase in success as the F = 1 constraint is relaxed. 
We will use this ratio to compare the behavior of several 
gates. We note that the matrices U corresponding to the 
curve in Fig. [2] remain unitary along the entire length 
of the curve and maintain the "Knill Form." The "Knill 
Form" is a empirical statement about the structure of the 
U matrix for optimal success with perfect fidelity in the 
case of the controlled-unitary and Toffoli gates, where one 
mode in each qubit maybe chosen to be non-interacting 
(up to local rotations) jl|, [l2|, [l3| . 

We also note that as mentioned above, this curve ob- 
tained by continuously deforming the F = 1 optimal so- 
lution is not guaranteed to represent the global maximum 
of the success rate for all values of F. Indeed it is not 
surprising that for large deviations from perfect fidelity, 
5 = 1 — F > .25, we may obtain solutions with higher 
success rate S for the same fidelity F, corresponding to 
other families of optimal solutions. Fig. [3] shows this for 
the case of the CNOT gate, in which we find two other 
families of solutions. One of these has success rates sig- 
nificantly higher than the "Knill Form" solution when 
S > .25 (F < 0.75), and the other has a fidelity-success 
curve that crosses that of the "Knill Form" family. While 
these results are interesting, the large S values corre- 
sponding to fidelities F < 0.75 make these alternative 
families of solutions less useful in practice as gates. 

As discussed in the previous section, the success rates 
shown in Fig. [3] and in all following figures are given by 



FIG. 3: (Color Online) The numerical results of fidelity- 
success optimization for the CNOT gate, extended to a larger 
domain of 5 = 1—F than in Fig. [2j and showing two additional 
one-parameter families of solutions that cross or surpass the 
"Knill Family" of solutions at finite 5. 



Success Rate of the CS(90°) vs. 8 

and the success rate divided by the two smallest singular values (NSV) 




FIG. 4: (Color Online) Black circles show numerical results 
for the success rate S as a function of S = 1 — F for the 
cs(90°) gate, obtained by continuously deforming the optimal 
F = 1 solution. The solid curve (Eq. ©) is a best fit to the 
data. The green diamonds and red squares show the success 
rate divided by the normalized smallest and second smallest 
singular values, respectively. 

The same logic that describes the shape of the fidelity- 
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success curve for cnot holds true for any LOQC gate, 
including cs(90°), the B gate, and the Toffoli gate. In 
Fig. i we show several results for cs(90°). The black 
circles show the success rates for cs(90°) found using our 
perturbative optimization method. The solid curve is the 
best fit to the theoretical form flJJ), and has the values 

S(S) = 0.05165 + .08368 5 1/2 + 0.23282 5 . (9) 

As in the case of cnot, the S = point confirms results 
found previously [H, [HI. The ratio Si/ So = 1.62 for the 
cs(90°) gate is larger than that of the cnot gate. 

In contrast with the cnot solution, where the opti- 
mal linear transformations U associated with the Knill 
family of solutions are all unitary, the optimal (N c + 
N a ) x (N c + N a ) = 6x6 matrix for the cs(90°) gate 
at F = 1 (5 = 0) has a four-fold degenerate maximum 
singular value and two smaller non-degenerate singular 
values [l2|. As discussed above, the maximum singu- 
lar value is always normalized to unity. This means in 
practice that the physical device implementing this op- 
timal solution requires two additional vacuum modes for 
unitary dilation [l|, 0, [l3j , for a total number of modes 
N c + N a + N v = 8. 

The behavior of the non-degenerate singular values 
when the perfect fidelity condition is relaxed can be seen 
in Fig. HJ As noted earlier, the black circles show the 
maximal success rate at each value of the fidelity, which 
may also be thought of as the maximal success rate di- 
vided by the largest normalized singular value (NSV), 
which is always normalized to unity. The green diamonds 
show the success rate divided by the smallest NSV at 
each point, and similarly, the red squares show the suc- 
cess rate divided by the second smallest NSV. Here we 
observe interesting behavior: as we move away from the 
perfect fidelity solutions by increasing 5, the singular val- 
ues converge. Now the second smallest NSV starts out 
within a few percent of unity, so it is difficult to deter- 
mine exactly where it merges with the maximal NSV. 
The data associated with the smallest NSV is more illus- 
trative: the green diamonds show that the smallest NSV 
merges with the other singular values at S ~ 0.03. At 
this point the optimal solution becomes unitary, and it 
remains unitary for larger 5, as indicated by the three 
data sets overlapping. This means that not only is our 
method optimizing the success rate at a given fidelity; 
it is actually optimizing the system size as well, as uni- 
tary dilation is no longer needed beyond the transition 
point. This reduction in system size and singular value 
convergence will be seen again in other gates. 

We remark that while the perturbative expansion (0 
works very well up until the NSV coalescence point at S « 
0.03, it clearly fails at larger values of 5. This is consistent 
with the fact that the structure of the optimal matrix U is 
different on the two sides of the transition. The solutions 
with different degeneracy of the singular value spectrum 
should be regarded as constituting separate families of 
optimal transformations, although the optimal success 
rate behaves continuously at the coalescence point, where 



a transition between the two families occurs. 

The final two-qubit gate we examine is the B gate 
found by Zhang et al. [14[, which in the Cartan decompo- 
sition has coordinates (tt/2, 7r/4, 0). The B gate is a max- 
imally entangling gate that has the unique property of 
only requiring two copies of the gate, along with generic 
single-qubit rotations, to produce an arbitrary two-qubit 
(5E/(4)) transformation. Any other two-qubit gate, in- 
cluding CNOT, requires three or more instances of itself 
to produce a general SU(4:) transformation in the stan- 
dard circuit model p], [l4| . Therefore the B gate can be 
thought of as the most general of all two-qubit gates. Pre- 
viously the authors found numerically that the B gate has 
a maximum success rate of S ~ 0.00717 at 5 = 0. Unlike 
the controlled-phase gates (equivalent to the controlled- 
unitary C X \J gates), which require only two single-photon 
ancillas, the B gate or any other two-qubit gate requires 
a third single-photon ancilla to achieve perfect fidelity. 
This makes the U matrix 7x7 instead of 6 x 6, while 
the A matrix remains the same size as the computational 
space of two qubits is unchanged. In Fig. [5] we see the 
results. 



Success Rate of the B-gate vs. 8 
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FIG. 5: (Color Online) Black circles show the maximal success 
rate as a function of 5 = 1 — F for the B gate, obtained by 
continuously deforming the optimal F = 1 solution. The solid 
curve (Eq. ([TO))) is a best fit to the data. The blue triangles, 
green diamonds, and red squares show the success rate divided 
by each of the three smallest singular values respectively. 

The B gate optimal solution at perfect fidelity is non- 
unitary, with three single-photon ancilla modes and three 
vacuum modes. As in the case of the cs(90°) gate, the 
singular values converge as S increases, with the first NSV 
merging with unity at S « 0.004 and the second merging 
at 5 « 0.055. The success rate for small S is given by 

S(S) = 0.0071 + 0.0308 S 1/2 + 0.0129 S , (10) 

making the key ratio S±/ Sq = 4.34 significantly larger 
than for cnot or cs(90°). Because the B gate solution is 
not in "Knill Form," it has full access to the 7x7 space, 
ignoring the vacuum modes, of the U matrix. This space 
is larger than that of either cs(90°) or cnot, where the 
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6 x 6 U matrix is constrained by the "Knill Form" to 
have two non-interacting modes near (5 = 0, giving the 
matrix an effective size of 4 x 4. We next examine the 
three-qubit Toffoli gate and compare the ratio Si / So for 
the Toffoli gate with the two-qubit gate results. 



IV. THE TOFFOLI GATE 

The Toffoli gate is the logical extension of the cnot 
gate to three qubits. It is sometime referred to 
as Controlled-Controlled-NOT, since the third qubit is 
flipped if and only if the first two qubits are both in the 
"on" state. The authors [13] found that the maximum 
success rate of the Toffoli gate at 5 = is S ~ 0.0034. 
Again, the success rate for imperfect fidelity should be- 
have as in Eq. ([7]) . The results of the perturbative relax- 
ation can be seen in Fig. [6] 



Success Rate of the Toffoli Gate vs. 8 

and the success rate divided by the smallest normalised singular values (NS V) 
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to unity at S « 0.075 and 5 « 0.15, respectively. This 
means that the Toffoli gate requires a total of 12 modes 
for its implementation at F = 1 (N c = 6 computational 
modes for the qubits, N a = 3 ancilla modes, and N v = 3 
vacuum modes), but we can reduce the system size to 
11x11 at S « 0.075 and further down to 10 x 10 at 
5 ~ 0.15. Just as in these case of the cs(90°) and B 
gates, our method is optimizing not just the success rate 
but also the physical size of the system. 

Fig.[3shows the success rate and curve fit as well as the 
behavior of the remaining singular values more clearly by 
dropping the data associated with the smallest singular 
value. 



Success Rate of the Toffoli Gate vs. 8 

and the success rate divided by the smallest normalised singular values (NSV) 
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FIG. 7: (Color Online) The results are as in Fig. [6] except 
that here we drop the data showing the smallest normalized 
singular value (NSV1) to provide better detail of the behavior 
of the others. 



FIG. 6: (Color Online) Black circles show the maximal success 
rate as a function of 5 — 1 — F for the three-qubit Toffoli 
gate, obtained by continuously deforming the optimal F — 1 
solution. The solid curve (Eq. (|lip ) is a best fit to the data. 
The blue triangles, green diamonds, and red squares show 
the success rate divided by each of the three smallest singular 
values respectively. 

As with the two-qubit gates shown above, we show the 
success rate as a function of fidelity using black circles, 
along with a best fit perturbative expansion for small 5. 
The Toffoli gate can be implemented with three single- 
photon ancilla modes. Since a three-qubit gate requires 
six computational modes, the linear transformation U is 
a matrix of size at least 9x9. Furthermore, the opti- 
mal solution at perfect fidelity has three non-degenerate 
singular values P, [l3j |. thus requiring 3 vacuum modes. 
It is also a "Knill Form" solution, in that one mode of 
each computational qubit is non-interacting, and just as 
in the case of cnot this form holds for the entire family 
of optimal solutions shown in the figure. Fig.[6]also shows 
the success rate divided by the three non-degenerate nor- 
malized singular values. Here we see that the second 
and third smallest normalized singular values converge 



The Toffoli gate success rate for imperfect fidelity be- 
haves as 

S(S) = 0.0039 + 0.0089 S 1/2 + 0.0258 S . (11) 

The ratio Si/ So has the value 2.28; this is squarely be- 
tween the values of 4.34 and 1.62 obtained for the B gate 
and the cs(90°) gate, respectively. This is not surprising, 
for while the Toffoli gate is a three-qubit gate, its solution 
obeys the "Knill Form," and the U matrix is effectively 
reduced from 9 x 9 to 6 x 6, ignoring the variable vacuum 
modes. This puts the results in line with our findings 
for the two-qubit gates. A summary of the results for 
all four gates appears in Table HI As seen in the table, 
gates requiring a larger number of interacting modes in 
the optimal solution at perfect fidelity have greatest po- 
tential for improved success rates as the F = 1 condition 
is relaxed. 

As with the cs(90°) and B gates, the analytic curve 
given by Eq. (O is valid in the region between 5 = and 
the point at which the largest non-degenerate singular 
value converges to the matrix norm. In the case of the 
Toffoli gate, this is the point S ~ 0.075. Beyond this 
point, the success rate jumps slightly above the curve, 
as compared with the cs(90°) gate, where the success 
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Gate 


Effective 

N c + N a 


Effective 

N c + N a + N v 


Si/ So 


CNOT 


4 


4 


1.03 


cs(90°) 


4 


6 


1.62 


Toffoli 


6 


9 


2.28 


B 


7 


10 


4.34 



TABLE I: The table summarizes results for four gates investi- 
gated in this paper. Here Effective N c -\-N a is the total number 
of interacting computational and ancilla photon modes in the 
optimal solution at perfect fidelity. For solutions obeying the 
Knill Form pi [I2L , N c equals the number of qubits; other- 
wise N c is twice the number of qubits. Effective N c + N a + N v 
adds the number of required vacuum modes in the optimal so- 
lution at perfect fidelity. Si /So (Eq. JB}) measures the trade 
off between fidelity and success rate as the perfect fidelity 
condition is relaxed: larger values of this ratio indicate larger 
relative increase in the success rate as fidelity is reduced. 



rate falls below the curve immediately after crossing the 
NSV coalescence point. We also note that the second 
largest non-degenerate NSV accelerates its approach to- 
wards unity once the largest non-degenerate NSV has 
converged, i.e., once the system size has been reduced 
from 12 x 12 to 11x11. This latter behavior is also man- 
ifested in the case of the B gate, although it is easier to 
observe in the case of the Toffoli gate in Fig. [71 

Finally, we observe numerical noise in the data, associ- 
ated with the optimization returning to the ridge of opti- 
mal solutions after a random perturbation at each step. 
The robustness of the results to perturbation strength- 
ens our claim that the fidelity-success curves shown in 
the figures, at least near S = 0, indeed correspond to the 



globally optimal solutions. 



V. CONCLUSION 

We have expanded on our previous knowledge base for 
optimization of LOQC gates. We have used new meth- 
ods to test the relatively unexamined area of LOQC gates 
with imperfect fidelity for both two-qubit and three-qubit 
gates. We have calculated the ratio Si/ Sq for four differ- 
ent gates and found that it is a useful parameter in char- 
acterizing the behavior of imperfect fidelity gates. Our 
results imply that gates involving many photon modes, 
which tend to have smaller success rates, benefit most 
from relaxing the prefect fidelity constraint. We have also 
shown that LOQC gates that are non-unitary at perfect 
fidelity (i.e., gates requiring additional vacuum modes for 
their implementation) need fewer and fewer resources as 
the fidelity is reduced, until the optimal solutions even- 
tually become unitary. Future directions include iden- 
tifying physical effects that can lead to these imperfect 
gates, calculation of the success rates for individual in- 
put states, such as the Bell states, and investigating the 
behavior of single quantum state generation. 
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